Diffractive Hadroproduction of Dijets and W's at the Tevatron Collider 
and the Pomeron Structure Function 
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Results from a phenomenological analysis of dijet and W hard diffractive hadroproduction at 
Fermilab Tevatron energies are reported. The theoretical framework employed here is a modi- 
fied version of the Ingelman-Schlein approach which includes DGLAP-evolved structure functions. 
Different from what has been achieved by the DESY ep HERA reactions, a reasonable overall de- 
scription of such diffractive hadron processes is obtained only when a complex, quark-rich Pomeron 
structure function is employed in the calculation. 
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I. INTRODUCTION 

Regge phenomenology is well known for providing a 
suitable and economical theoretical framework for the 
description of soft hadron diffractive processes at high 
energies. Single Pomeron exchange plus a few secondary 
Reggeon contributions are enough to describe a variety 
of hadronic reactions (see, for instance, and references 
therein) . 

The situation becomes much more intricate when one 
considers hard diffractive processes by which, according 
to the Ingelman-Schlein (IS) picture Q, the Pomeron 
structure itself is probed. Difficulties arise when one tries 
to obtain a unified description for diffractive processes 
starting with both electron- or positron-proton (ep) and 
antiproton-proton (pp) collisions, respectively, studied at 
the DESY ep HERA and Tevatron colliders. Although 
several theoretical approaches have successfuly been em- 
ployed to describe different aspects of hard diffraction 
revealed by the ep HERA reactions some of them 
based on Regge theory, diffractive hadroproduction con- 
tinues to be one of the most challenging topics in hadron 
dynamics. 

In effect, most of the theoretical approaches which are 
able to describe the HERA data are not readily translat- 
able to diffractive hadron physics. Models based on the 
Regge theory, in particular, are presumably affected by 
a lack of validity for QCD factorization in the hadronic 
diffraction domain Q. In spite of these difficulties, such 
models establish the phenomenological picture most of 
the event generators currently employed in data analy- 
sis of hard diffractive hadroproduction are based upon. 
Probably this is so because these models have been able 
to provide an effective description for such processes. In 
fact, that is one of the underlying assumptions of the 
present paper. 

This Brief Report is a sequel to a previous work |5j 
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in which we have tried to perform a global analysis by 
a modified version of the IS model, including processes 
initiated by both ep and pp collisions. By that time, the 
available data were not so stringent such that one could 
speculate about an unique model to be sufficient. Since 
then, more and more precise data of diffractive deep in- 
elastic scattering (DIS) have imposed severe restrictions 
on the Pomeron structure function making more evident 
the impossibility to readily transfer the partonic densities 
so obtained to hadronic process calculations. 

Our group's analysis of the diffractive DIS data have 
shown that, at low values of the QCD evolution scale, 
the Pomeron is predominantly composed by gluons with 
a hard distribution, in agreement with other studies (see 
[6( and references therein). This result is corroborated 
by the calculations of diffractive cross sections for photo- 
production 7] and electroproduction 8] of dijets, which 
yield good agreement with the data once a hard gluonic 
Pomeron is assumed. 



The advent of new data produced by the DO Collabo- 
ration have motivated us to perform a new analysis, 
this time restricted to the Tevatron data. 



Thus, we report here results of a study on diffrac- 
tive hadroproduction of dijets and W's by using the IS 
model, with Dokhshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution [lTJ included, but disconnected from 
the HERA data analysis just mentioned [(j. In fact, in 
spite of our efforts it was impossible to reconcile both 
analyses. We take this failure to produce a unified de- 
scription as an additional indication of the theoretical 
problems alluded before. 

The interesting point, however, is that if one takes 
the Pomeron as predominantly composed by quarks at 
a low QCD evolution scale, a reasonable overall descrip- 
tion of diffractive hadroproduction data is achieved. This 
is what is shown below. 
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II. CROSS SECTIONS 



and 



Our starting point is the generic cross section for a pro- 
cess in which partons of two hadrons, A and B, interact 
to produce jets (or Ws), A + B —> Jets (W) + X, that 
is 

d(T = ^ fa/A(Xa,^ 2 ) dx a fb/B^b,^) dx b 
a,b,c,d 

da„ 



>ab^cd(W) , r 

x ^ — - at. 

dt 



(1) 



From this very basic expression we derive all of the oth- 
ers necessary to describe the specific processes we are 
interested in. 



A. Diffractive dijet production from a single 
Pomeron exchange 

In the case of dijet production, the cross section can be 
put in terms of the dijet rapidities (77, rj') and transversal 
energy Et- 



dr)J a ♦ JE Tm ,„ T in' ' 

' ' J J partons ^Tmxn '/ m in 



x x a f a / A (x a ,fi ) x b f b/B (x b ,fj, ) I —j (2) 



d 4\ 



dt J 



where 



x a = ^( e -" + e -"') ) Xb = ^(g-J + gi'), (3) 
V s V s 



with 



In- 



< 77 < In 



<s- E T e-v 



^fs-E T e~2 



(4) 
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T max 



(•5) 



being that Et min an d the 77 range are determined by the 
experimental cuts. 

Equations. ©-© express the usual leading-order 
QCD procedure to obtain the non- diffractive dijet cross 
section (next-to-leading-order contributions are not es- 
sential for the present purposes; see Ref. In or- 
der to obtain the corresponding expression for diffrac- 
tive processes, we assume that one of the hadrons, say 
hadron A, emits a Pomeron whose partons interact with 
partons of the hadron B. Thus the parton distribu- 
tion x a f a /j\(x a ,^, 2 ) in Eq. is replaced by the con- 
volution between a putative distribution of partons in 
the Pomeron, (3f a /w(P, M 2 )) and the "emission rate" of 
Pomerons by A, /jp(xp,i). This last quantity, fjp(xjp,t), 
is the so-called Pomeron flux factor whose explicit for- 
mulation in terms of Regge theory is given ahead. The 
whole procedure implies that 

Xafa/A{Xa,^) = Jdx w J d(3 J dt f w (x w ,t) 

x /3/a/i>(/3,M 2 H(/3-— ), (6) 

Xjp 

and, defining g(xsp) = J_ dt fw(xw,t), one obtains 



Xafa/A(Xa, V 2 ) = I dx w g(x W ) — f a/w ( — , (I 2 ) . (7) 
J Xjp Xjp 

By inserting the above structure function into Eq. J5J one 
obtains the cross section for diffractive hadroproduction 
of dijets via a single Pomeron exchange as 



d(TSPE \ _ yi I I d l I dx W g(x w ) Pafa/TpiPa,^ 2 ) X b f b /p{x b ,fi 2 ) (^Z^p±\ (g) 



where (3 a — x a /xsp with x a and Xb given by Eq. 0, and 
x TP mill an d x w max established by the experimental cuts. 

B. Diffractive Hadroproduction of W 

diffractive production is here considered by the re- 
action p + p — > p + W{— * e v) + X. It is assumed that a 
Pomeron emitted by a proton in the positive z direction 



interacts with a p producing W that subsequently de- 
cays into e v. The detection of this reaction is triggered 
by the lepton (e + or e~) that appears boosted towards 
negative rj (rapidity) in coincidence with a rapidity gap 
in the right hemisphere. 

By using the same concept of the convoluted structure 
function, the diffractive cross section for the inclusive 
lepton production for this process becomes 



3 



dasp \ 
drje J 



lepton a ft 



^2 J g ( Xw ^ J dEr /i/ff(a;a,M 2 ) fb/p{x b ,y?) 



6 s r 



w 



(9) 



r 



where 



M w e* 1 ' 
(y/s Ip) 



M w e-v* 



A±V(A 2 -1) 



and 



t = — Ex M w 



A + ^{A 2 -l) 



(10) 
(11) 
(12) 



with A = Mw/^Et- The upper signs in Eqs. (|TU1) and 
(|ll|l refer to W + production (that is, e + detection). The 
corresponding cross section for W~ is obtained by using 
the lower signs and i <-» u (see the Appendix in Q)- 

C. The Pomeron Flux Factor 

An important element of this approach is the Pomeron 
flux factor, introduced in Eq. 10. It has some peculiar 
aspects that deserve to be pointed out. 

First of all, the expression for this term was originally 
proposed to be taken from the invariant cross section of 
(soft) diffractive dissociationprocesses as it is given by 
the Triple Pomeron model |2j. The rationale for that 
can be put in terms of an analogy with the photon flux 
factor, this one derived from QED. The basic idea is that, 
similarly to what happens to the electron (or positron) in 
photoproduction, the proton in a diffractive interaction 
is scattered at very small angles and practically does not 
take part in the effective reaction. Analogously to the 
emission of photons and to the idea of equivalent photon 
flux defined in QED, one can think of hadron diffraction 
in terms of Pomeron emission and the "Pomeron flux 
factor" . This picture (and the IS model as a realization of 
it) has been successfully employed to the hadron vertex in 
some HERA diffractive processes, such as leading baryon 
production and diffractive DIS |6j], photoproduction [7j, 
and electroproduction [8|. 

However, such an approach is affected by a problem 
which is mostly concerned with its energy dependence. 
As it is theoretically well known from very long, the 
Triple Pomeron integrated cross section violates unitar- 
ity 0], although its and t dependences seem to be 
in good agreement with the available data I n order 
to overcome this unitarity violation issue, we follow here 
the "renormalization" procedure originally proposed in 
Hjl and further discussed in fiji, that is 



f{xw,t) 



4T~ It=of( x ^^) dx v dt 



(13) 



For the "unnormalized" flux factor f(xw, i), we take 
the Donnachie-Landshoff parametrization 15], 



where Fi{t) is the Dirac form factor, 



[Am 2 - 2.79t) 1 
(4m 2 - t) (1-^)2 



(14) 



(15) 



Notice that, by choosing the renormalization proce- 
dure, Po does need to be specified since it is crossed out as 
well as the other constant factors appearing in Eq. I|14|l . 
Yet about this equation, our choice for the Pomeron tra- 
jectory has been a%>(t) = 1.2 + 0.25 t, which is compatible 
with both Tevatron and HERA data. 



D. The Pomeron Structure Function 

The Pomeron structure function has been established 
as a three-flavor quark singlet at the initial scale, chosen 
to be Qq = 2 GeV 2 , with the gluon component being gen- 
erated by DGLAP evolution. Thus, no initial gluon dis- 
tribution has been assumed. The parametrization used 
for the initial quark distribuition was 

/?£(/?, Ql) = [At exp(-A 2 /3 2 ) +B 1 {1- f3) B *} (3 - 001 
+ Ci exp[-C 2 (l - /3) 2 ](1 - /3) - 001 , (16) 

which includes different amounts of soft, hard, and su- 
perhard profiles according to the chosen parameters. The 
results presented below were obtained with the following 
parameters: A\ — 4.75 and A2 — 228.4 for the soft part, 
B\ = 1.14 and B2 — 0.55 for the hard one, and finally 
C\ = 2.87 and C2 = 100 for the superhard term. 

Wherever necessary, DGLAP evolution of the Pomeron 
parton densities has been processed by using the program 
QCDNUM 0. For the proton (or antiproton, when 
was the case), the parton densities were taken from the 
parametrizations given in Ref. [Tt) . 



III. RESULTS AND DISCUSSION 

In the following, we present our predictions for hard 
diffractive production of W's and dijets based on the pre- 
vious discussion. These predictions are compared with 
experimental data from Refs. P. IT8llT^ | in Tables I-III. 

In Table I, the difficulty in obtaining a perfect and 
simultaneous description of both W and dijet production 
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TABLE I: Data versus model predictions. Diffractive W's 
and dijets were measured at y*s = 1800 GeV by the CDF 
Collaboration [Til Il9|. In both cases, i r < 0.1 and Er min = 
20 GeV. For the case of W production, Ex min refers to the 
detected lepton while for dijet production it refers to the de- 
tected jet. 



Yield 


Rapidity 


Data (%) 


Model 


W 


-1.1 <n e < 1.1 


1.15 ±0.55 


0.35 


jj 


-3.5 < rjj < -1.8 


0.75 ±0.10 


0.72 



ratios are calculated: ratios between rates at different 
energies but at the same rapidity range and the reverse, 
ratios between rates taken at the same energy but differ- 
ent rapidity ranges. The better agreement here could be 
attributed to the fact that these ratios would cancel the 
normalization and acceptance effects to some extent. 

In summary, we have shown that it is possible to ob- 
tain a reasonable overall description of hard diffractive 
hadroproduction by a model based on the Ingelman- 
Schlein approach once a quark-rich Pomeron structure 



TABLE II: Data versus model results corresponding to the 
DO experiment. The experimental data are from Ref. [jj and 
the model calculations were performed with Er min = 15 GeV 
for y/s = 1800 GeV and E Tmm = 12 GeV for -Js = 630 GeV. 
In both cases, x w < 0.1. 



(GeV) 


Rapidity 


Data (%) 


Model 


1800 


\V\ > 1-6 


0.65 ± 0.04 


0.90 


1800 


M < i.o 


0.22 ±0.05 


0.37 


630 


\V\ > 1-6 


1.19 ±0.08 


1.80 


630 


M < i.o 


0.90 ± 0.06 


0.98 



is evident. The situation is much better, however, when 
one considers only jets. Besides the agreement exhibited 
in Table I for the CDF experiment, consistency is also 
found with the DO results (Tables II and III). 

In Table II both forward and central dijet production 
at two energies are considered. For all cases, one sees that 
the model predictions are close to the data, but slightly 
above. This sort of discrepancy is expected since effects 
of experimental acceptance were not taken into account 
in these predicitions. Such effects would certainly reduce 
these theoretically predicted rates, but it is difficult to 
estimate to what extent. 

In Table III is where the agreement between theory 
and data is generally better. In this case, two kind of 



TABLE III: Experimental ratios versus model results corre- 
sponding to the DO experiment. Data are from Ref. [jj. 



Ratios 


Data (%) 


Model 


630/1800 


M>1.6 1.8 ±0.2 


2.0 


630/1800 


|r/| < 1.0 4.1 ±0.9 


2.7 


1800 GeV 


\r}\ > 1.6/M < 1.0 3.0 ±0.7 


2.4 


630 GeV | 


r?| > 1.6/|r?| < 1.0 1.3 ±0.1 


1.8 



function is assumed and its DGLAP evolution performed. 
This result, i.e. the predominance of quarks in the 
Pomeron "valence" distribution, already obtained in 
is in conflict with the parametrizations independently es- 
tablished from HERA data H Q, ||. This discrepancy 
may be seen as an additional indication of factorization 
breaking £| in hadronic diffraction. However, if that is 
the real reason, it is quite intriguing that the consistence 
between the data and theory shown here is possible at 
all. 
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